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Abstract 

We study loop operators of 5d N = 1 SYM in Q background. For the case of U(l) 
theory, the generating function of correlation functions of the loop operators reproduces 
the partition function of melting crystal model with external potential. We argue the 
common integrable structure of 5d Af = 1 SYM and melting crystal model. 
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1 Introduction 



It is shown in pQ that the Seiberg-Witten solutions [2] of 4c? Af = 2 supersymmetric gauge 
theories emerge through random partition, where Nekrasov's formulas [U [3] for these gauge 
theories are understood as the partition functions of random partition. The integrable structure 
of random partition is elucidated in [I], and thereby the integrability of correlation functions 
among single-traced chiral observables is explained. Such an extension of the Seiberg-Witten 
geometries also becomes attractive to understand 4d Af = 1 supersymmetric gauge theories by 
providing a powerful tool [5]. 

Integrable structure of melting crystal model with external potential is clarified in [6] . Melt- 
ing crystal model, known as random plane partition has a significant relation with 5d A" = 1 
supersymmetric gauge theories. Nekrasov's formula for these gauge theories can be retrieved 
from the partition function of melting crystal model [7], where the model is interpreted as a 
g-deformed random partition. It is argued [6] a relation between loop operators of 5d Af = 1 
supersymmetric Yang-Mills (SYM) and external potentials of the melting crystal model. 

We start Section 2 with providing a brief review about 5d A" = 1 SYM in Q background 
[8]. We introduce loop operators of this theory. Computation of correlation functions among 
these operators is discussed. Generating function of the correlation functions of U(l) theory 
reproduces the partition function of the aforementioned melting crystal model. In Section 3 
we discuss a common integrable structure of 5d Af — 1 SYM in Q background and melting 
crystal model for the case of the U(l) theory. In Section 4 we present an extension of the 
Seiberg-Witten geometry of the U(l) theory by using the loop operators. 

2 Loop operators of 5d J\f = 1 SYM in Q background 

We first consider an ordinary 5d Af = 1 SYM on 1R 4 x S 1 . Let E be the SU{N) -bundle on M 4 
with 02(E) = n > 0. A gauge bundle of this theory is the S*?7(A^)-bundle tt*E on R 4 x S l pulled 
back from IR 4 . tt is the projection from M 4 x S 1 to 1R 4 . All the fields in the vector multiplet are 
set to be periodic along S 1 . The bosonic ingredients are a 5d gauge potential Am(x, t)dx M and 
a scalar field ip(x,t) taking the value in su(N). These describe a 5d Yang-Mills-Higgs system. 
The gauge potential can be separated into two parts A^x, and A t (x, t)dt, respectively the 
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components of the M 4 - and the S^-directions. Let Ae be the infinite dimensional affine space 
consisting of all the gauge potentials on E. A^(x,i)dx^ describes a loop A(t) in Ae, where 
the loop is parametrized by the fifth-dimensional circle. As for A t (x,t), together with (p(x,t), 
the combination A t + i<p describes a loop <p(t) in S1°(M. 4 , adE ® C), the space of all the sections 
of adE <S> C, where adE" is the adjoint bundle on M 4 with fibre su(N). Taking account of the 
periodicity, the same argument is also applicable to the gauginos. The vector multiplet thereby 
describes a loop in the configuration space of the 4d theory. In the case of the Yang-Mills-Higgs 
system, the loop A(t) gives a family of covariant differentials on E as d^A — d + A(t). For the 
loop (pit), since it involves A t (x,t), it becomes convenient to introduce the differential operator 

W (f) = ^ + 0(t). (2.1) 
2.1 5dN=l SYM in Q background 

Via the standard dimensional reductions, 6d M — 1 SYM gives lower dimensional Yang-Mills 
theories with 8 supercharges, including the above theory. Furthermore, the dimensional re- 
ductions in the Q background provide powerful tools to understand these theories [8]. The Q 
background is a 6d gravitational background on IR 4 x T 2 described by a metric of the form: 
ds 2 = X]J = i(^' 1_ Efl=5 6 ^adx a ) 2 + ^ a=5 6 {dx a ) 2 , where two vectors V^, V§ generate rotations 
on two-planes (x 1 , x 2 ) and (x 3 , x 4 ) in M 4 . By letting V± = x 2 -£^ — x 1 -^ and V2 = x 4 ^ — x 3 -^, 
they are respectively the real part and the imaginary part of the combination 

V ei , e2 = e 1 V 1 + e 2 V 2 , ei,e 2 GC (2.2) 

The above combination is expressed in component as V €l>e2 = £l fJi v x u . 

To see the dimensional reduction in the ^-background, we first consider the bosonic part 
of the 5d SYM. The corresponding Yang-Mills-Higgs system is modified from the previous one. 
However, the system is eventually controlled by replacing Ti{t) with 

n eue2 (t)=H(t) + }C eue2 (t). (2.3) 

Here /C ei>£2 (t) is an another differential operator of the form [9] 

Ke^{t) = v^ t2 d mil + larjr , (2.4) 
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where denote the 5*0(4) Lorentz generators of the system. This operator generates a 
T 2 -action by taking the commutators with d A (t) an d 7~L{t)- For instance, we have 

[d A( t),lC eue2 (t)] = -iv ei , e2 F A( t) . (2.5) 

The right hand side is precisely the transformation brought about on Ae by the infinitesimal 
rotation 5x^ = —V£ . 

The supercharges Q aa and Q% are realized in a way different from the case of e± = e 2 = 0. 
Note that we use the 4c? notation such that a, a and a denote the indices of the Lorentz 
group SU(2)l x SU(2)r and the R-symmetry SU(2)j. By the standard argument, we may 
interpret the 5d SYM as a topological field theory. Actually, by regarding the diagonal SU (2) 
of SU{2) R x SU{2) I as a new SU(2) R , we can extract a supercharge that behaves as a scalar 
under the new Lorentz symmetry. We write the scalar supercharge as Q ei ,e 2 - The gaugino 
acquires a natural interpretation as differential forms, r)(x, t), ip^x, t) and ^ u (x,t). These give 
fermionic loops, r](t),ip(t) and £(£). The main part of the Q-transformation takes the forms 

Q £1 , £2 A(t) = m , Q ei , e M) = [d A{ t),n ei , e2 (t)} , (2.6) 
Q £l)£2 W £l)£2 (t) = 0, (2.7) 

where ip(t) is a fermionic loop in f2 1 (M 4 , &dE). 

2.2 Loop operators and their correlation functions 

Taking account of the relation (p(x,t) = A t (x,t) + i(p(x,t), the following path-ordered integral 
provides an analogue of a holonomy of the gauge potential. 

W^ix-MM) = Pe-^ 1 ^') , (2.8) 

where the symbol means the path-ordered integration, more precisely, it is defined by the 
differential equation 

(-L + <i>(x,t 1 ))W<- \x-,t 1 ,h) = 0, WM(x;t 2 ,t 2 ) = 1. (2.9) 
at 1 

The trace of the holonomy along the circle defines a loop operator as 

{0) {x) =TyW (0 \x;R,0), (2.10) 
3 



where R is the circumference of S . The above operator is an analogue of the Wilson loop 
along the circle. Unlike the case of e± = e 2 = 0, it is not Q-closed except at x — 0. To see this, 
note that the Q-transformations (12.61) and (12. 7p imply Q eie2 (j)(t) = — ip(t). By using this, 
we find 

Q^ t2 O i0 \x) = j* dtxTi^W^fa RM)iv ei ^{xM)W^{x- t x ,0)}. (2.11) 

Since the right hand side of the above formula vanishes only at x — 0, this means that O^ \x) 
becomes Q-closed only at x — 0. 

The above property may be explained in terms of the equivariant de Rham theory. To see 
this, let us first generalize the path-ordered integral (12.81) by exponentiating the combination 
Fa{£) — i>{t) + (fi(t) in place of 0(t) as 

W{x- h,t 2 ) = p e -/4 ld *( F ^)-^)+^))w , (2.12) 

where the right hand side is given by a differential equation similar to (12. 9p . This means that 

W has the components, according to degrees of differential forms on IR 4 , as 

..._(_ jy( 4 ) ; where the indices denote the degrees. We generalize the loop operator ( I2.10p as 

0{x) = TtW(x;R,0) . (2.13) 

This also has components as O = + + ■ • • + O^. Eq. ( 12.111) can be now expressed as 
Q ei , t2 0^ = Lv Cl , C2 0^\ This is actually the first equation among a series of the equations that 
obey. Such equations eventually show up by expanding the identity [10] 

(d eii6a + Q ei , 6a )O(ar) = 0, (2.14) 

where <i £li£2 = d — Ly ei e2 is the T 2 -equivariant differential on IR 4 . 

We can also consider the loop operators encircling the circle many times. Correspondingly 
we introduce 

O k (x) = TrW(x;kR,0), k = l,2,--- (2.15) 

These satisfy 

{d euea +Qe uea )O k (x) = 0. (2.16) 
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Let us examine the correlation functions ( Y[ a J K 4 @k a ) €1 ' £2 - Since the integral / R4 = 
Jr4 i s Q-closed by virtue of the formula (|2,16p . these can be computed by a supersymmetric 
quantum mechanics (SQM) which is substantially equivalent to the 5d SYM as the topological 
field theory. Such a SQM turns to be T 2 -equivariant SQM on Ai n [3], where M. n is the 
moduli space of the framed n instantons. The Q-transformation (12. 6p is converted to the 
supersymmetry of the quantum mechanics 

Qe^eMt) = X(t) , Qe 1 ,e 2 X(t) = + V ei , ea (m(t)) , (2.17) 

where V tltt2 is the Killing vector induced by the variation 5 A = iv ei , e on Ae- The combi- 
nation Fj^m — ip{t) + (j)(t) can be identified with a loop space analogue of the T 2 -equivariant 
curvature T ei ^ 2 of the universal connection, where the universal bundle becomes equivariant 
by the T 2 -action on A E x E 4 . 

In the computation of the correlation function, by virtue of the supersymmetry (I2.17p . 
only the constant modes mo, Xo contribute to the observable, and the above combination pre- 
cisely becomes T ex ^ 2 [8]. This means that Ok(x) truncates to the equivariant Chern character 
Tr e - fcR ^i,e 2 . Thus we obtain the finite dimensional integral representation 

( II / T'" = [ A T >{Rt ex>£2 , M n ) J] / Tre-^^- . (2.18) 

X a Jm. 4 1 n-mstanton (2lliR) 2 J Mn a 

where A T 2(- ,A4 n ) is the T 2 -equivariant A-genus of the tangent bundle of A4 n , and t ei e2 is a 
generator of T 2 that gives the Killing vector V eii<E2 - 

Introducing the coupling constants t = (ti, t2, ■ ■ ■ ), the generating function of the correlation 
functions is given . Since ra-instanton contributes with the weight 

(RA) 2nN , where A is the dynamical scale, letting Q = (RA) 2 , we can express the generating 
function as 

Z eut2 (t) = TQ nN (e^^°A eUe2 . (2.19) 

' * \ / n—instanton 

n=0 

2.3 Application of localization technique 

The right hand side of the formula (12.181) is eventually replaced with a statistical sum over 
partitions. To see their appearance, note that the integration localizes to the fixed points of 
the T 2 -action. However, the fixed points in M. n are small instanton singularities since the 
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variation 5 A = —lv £i , 2 Fa vanishes there. These can be resolved by instantons on a non- 
commutative M 4 . Applying such a regularization via the non-commutativity, the fixed points 
get isolated, so that they are eventually labelled by using partitions [TTj . 

A partition A = (Ai, A2, • • •) is a sequence of nonnegative integers satisfying A; > A*+i for 
all z > 1. Partitions are identified with the Young diagrams in the standard manner. The size 
is defined by |A| = J2i>i which is the total number of boxes of the diagram. 

Let us describe the formula (12. 18ft for the U(l) theory. The relevant computation of the 
localization can be found in [HI [12]. We truncate ei^ as —e\ = €2 = ih, where h is a positive 
real parameter. Consequently, the formula becomes a g-series, where q = e~ . The fixed 
points in Ai n are labelled by partitions of n. The equivariant A-genus takes the following form 
at the partition A of n: 

(2mR)- 2n A T >(Rt^ lh M n ) a = H n (^ ) (Y[Ks)) q^Sxtf)*, (2.20) 

where h(s) denotes the hook length of the box s of the Young diagram A, and s\(q p ) is the 
Schur function s\(xi, X2, ■ • • ) specialized to Xi = q l ~^. Similarly, the fixed points in M. n x M 4 
are labelled by partitions of n. Denoting them as (A, 0), the equivariant Chern character takes 
the form Tr e ~ kR:F -ih,ih | = fc (A), where Ok(X) is given by 

00 

O k (X) = (1 - q- k ) Y^{q k(X ^ l+1) - g fc( " i+1) } + 1 . (2.21) 

i=l 

The above functions have been exploited in [H [13] from the Ad gauge theory viewpoint. By 
taking account of (12.201) and (12.211) . the formula (12.181) becomes eventually as 

n / y lh,lh = (-r e i^^ p ) 2 n^o ka (x) . (2.22) 

/ n—instanton 

a l-M= n a 

Although we have not taken into account, the Chern-Simon term can be added to a 5d 
gauge theory, with the coupling constant being quantized, in particular, for the U(l) theory, 
m = 0, ±1. It modifies the right hand side of ( 12.221) by giving a contribution of the form 
(— ) m l A lg a , for each A [7j. Hereafter, we consider the case of the U(l) theory having the 
Chern-Simon coupling, m = 1. The corresponding generating function becomes 

= J2Q lXl s x (qn 2 e h - 2 ^--^ W . (2.23) 

A 
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3 Integrability of 5d Af = 1 SYM in Q background 



We can view the generating function (12.231) as a g-deformed random partition. To see this, 
note that the 4<i limit R — ** makes q = e~ Rh — > 1, the Boltzmann weight takes at this limit, 
the form (A//i) 2 ' A ' (H s eA > which is the standard weight of a random partition. It can be 

also viewed as a melting crystal model, known as random plane partition. The corresponding 
model is studied in [6] as a melting crystal model with external potential, where the Chern 
characters Ok correspond precisely to the external potentials. 

3.1 Melting crystal model 

A plane partition 7r is an array of non-negative integers 

Til 7Ti2 VT 13 • • • 
7T 2 l 7T 2 2 7T 2 3 • • • 

(3.1) 

7T31 7T32 ^33 

satisfying 7Ty > ir i+ ij and 7Ty > 7r^ +1 for all z,j > 1. Plane partitions are identified with the 
3d Young diagrams. The 3d diagram tc is a set of unit cubes such that 7i~jj cubes are stacked 
vertically on each -element of 7r. Diagonal slices of 7r become partitions, as depicted in 




(a) (b) 

Figure 1: The 3d Young diagram (a) and the corresponding sequence of partitions (b). 

Fig.l. Denote ir(m) the partition along the m-th diagonal slice, where m G Z. In particular, 
7r(0) = (71*11, ^221 • ■ ■ ) is the main diagonal one. This series of partitions satisfies the condition 

< tt(-2) -< tt(-1) -< tt(0) >- tt(1) >- ?r(2) > , (3.2) 
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where \l y v means the interlace relation; \i y v •<=>■ \x\ > V\ > \ii > v<i > A*3 > • • • • 

The hamiltonian picture emerges from the above relations, by viewing a plane partition as 
evolutions of partitions by the discrete time m. Eventually it is described [2] by using 2d free 
complex fermions if}, if}*. We may separate the relations (13.21) into two parts, each describing 
the evolutions for m < and m > 0. These two types of the evolutions are realized in the 2d 
CFT by using operators G± of the forms [l l 



G± = e^ =1 M^) ±fe , (3.3) 

where J±k = Yl^L-co '■ 4 , ±k-n4>n '■ are the modes of the 2/(1) current. 

Using the free fermion description, one can express the generating function as 

Z%nit) = <0|G+Q L W^^Gqo> , (3.4) 

where Lq = Xl^L-oo n '■ ^-n^Pn '■ is an element of the Virasoro algebra. The loop operators Ok 
are converted to operators Ok in the above representation. They are fermion bilinears given by 

Ok = (l-q- k ) : : + 1 > ( 3 - 5 ) 

n=— oo 

3.2 The integrable structure 

The fermion bilinears 0% can be regarded as a commutative sub-algebra of the quantum torus 
Lie algebra realized by the free fermions [6]. The adjoint actions of G± on the Lie algebra 
generate automorphisms of the algebra. Among them, taking advantage of the shift symmetry, 
the representation (13. 4p can be eventually reformulated [B] to 

Z% h (t) = (0| e ^£r=iR fe (i-<r fc )^* gdu{i) e ^Er=i(-)W)^-* | ) . (3.6) 
In the above formula, gl dU is an element of GL(oo) of the form 

where W = = ^^L.^ n 2 : if}- n if> n : is a special element of algebra. The loop operators 
Ok are converted to Jk or J_ fc in (13. 6p . These two are actually equivalent in the formula, since 
gl dU ^ satisfies [6] 

j k g^um = g 5du(i)j_ k ; for k > ( 3 8 ) 



Viewing the coupling constants t as a series of time variables, the right hand side of (I3.6P is 
the standard form of a tau function of 2-Toda hierarchy |15j . However, by virtue of (13.81) . the 
two-sided time evolutions of 2-Toda hierarchy degenerate to one-sided time evolutions. This 
precisely gives the reduction to 1-Toda hierarchy. Thus the generating function becomes a tau 
function of 1-Toda hierarchy. 

4 Extended Seiberg-Witten geometry of 5d theory 

We consider the field theory limit of the U(l) theory, which is achieved by letting H — > and 
amounts to the thermodynamic limit of the melting crystal model. The system is described by 
the prepotential J-^°'(t; A, R). From the integrable system viewpoint, may be interpreted 
as a dispersion-less tau function, since the generating function is substantially a tau function 
of 1-Toda hierarchy and J 7 ^ gives the leading order part of the h expansion of \ogZ^\^ ih {t) To 
obtain the semi-classical solution, one actually needs to solve the related variational problem, 
which is reformulated as a Riemann-Hilbert problem. This issue is treated in [TO] . 

4.1 Seiberg-Witten curve of U{\) theory 

Let us present the Seiberg-Witten curve for the U(l) theory. We first employ the following 
curve piU [17] : 

where (3 is a real parameter. Cp is a double cover of the cylinder C* = C/^p, with a cut / along 
the real axis on the Riemann sheet. The coupling constants t determine (5 as (3 — (3(t). To see 
this, let us introduce a meromorphic differential of the form 

= {l - -R 2 J2k 3 t k M k (z)}dlogy, (4.2) 

k=l 

where Mk(z) = Y^ n =o dk-n{(3)e~ nRz ■ The coefficients d n {(3) are given in the asymptotic expan- 
sion 

oo 

V^l - pe Rz ) 2 - (2RAe Rz ) 2 = J2 d n {P)e nRz , Mz -> -oo. (4.3) 

n=0 
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Finally, solving the Riemann-Hilbert problem, (3 is determined by the condition [10] 

<j) zdty = (C: a contour encircling J anticlockwise). (4-4) 

4.2 Vevs of the loop operators 

The vev of the loop operators Ok can be represented by using an analogue of the Seiberg-Witten 
differential. Eventually, the vev can be organized to the contour integral 

d^{t;A,R) r //n . -kR I _ kRzja , AK . 

y ' = hm(O k ) = —— <b e kRz dS, , 4.5 

dtk n-^o 2m J c 

where dS = S'(z)dz is an analogue of the Seiberg-Witten differential. S'(z) is given by the 

indefinite integral 



z 



S\z) = J . (4.6) 

The contour integral in the right hand side of H4.5[) can be converted to a residue integral. 
Actually, by using coordinate Z = e~ Rz , we obtain 
dF^(t;A,R) 
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